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In order to study the effect of tensor force, we plan to perform a (p,d) reaction with a 400-1600 MeV proton 
beam on the RIBLL2 at Lanzhou. Based on the experimental conditions of RIBLL2, a Monte Carlo method is 
used to simulate (p,d) reaction process. The distributions of primary beam and scattered deuterons are given 
on the target and at the F1, F2, F3, and F4 positions. Considering the yield of the deuteron, to separate the 
produced deuteron from the background particles, the target thickness is optimized. To obtain a clear particle 
identification spectrum by time of flight (TOF) and energy loss (AF), the distance between the two detectors, 
as well as the energy and timing resolution of detectors are simulated. As a result, the distance between F and 
F; is fit for that of TOF. After taking into account the particle distribution at F4, both sizes of Multi-Wire Drift 
Chamber, namely MWDC2 and MWDC3, are selected to be 50 mm x 50mm. 
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I. INTRODUCTION 


Tensor forces play an important role in the pion exchange 
interactions that provide the most important attraction in nu- 
clear forces [1-3]. The importance of the tensor forces has 
been demonstrated to reproduce the properties of nuclear mat- 
ter, as well as explain the binding energy [4] of deuteron and 
alpha particles. Admixture of D-wave in those nuclei pro- 
vides clear evidence of the tensor interactions [5]. Recently, 
theoretical calculations on °~!'Li that include explicitly the 
tensor interactions have pointed out [6] the importance of the 
tensor interactions in understanding the structure of those nu- 
clei to predict high momentum components in ground states. 

Experiments using the electron [7, 8] or proton-induced [9] 
knockout reaction have been performed to probe the tensor 
correlations in nuclei from !?C to ?°°Pb. However, due to 
other correlations such as short-range repulsion, it is dif- 
ficult to isolate the tensor effects unambiguously in these 
experiments. 

In fact, theoretical calculations [1, 3, 10] have predicted en- 
hanced momentum distributions at around 2fm~! due to the 
tensor interactions, mixing large orbital angular momentum 
states through the D-wave. The 2 fm~! has ~400 MeV/c mo- 
mentum transfer and corresponds to a 800 MeV proton beam. 

We aim to explore the effects of tensor force on 16O nu- 
cleus. In the shell model, the ground state of '°O consists 
of eight protons and eight neutrons filling up the 1s, /2, 1p3/2 
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and 11/2 orbits. Since the tensor interactions induce changes 
in the total orbital and spin angular momenta by |AL| = 2 
and |AS] = 2, the positive-parity states of '!°O can be reached 
through direct one-neutron pick-up reaction, only if the '°O 
ground state has an admixture of 1d5/2 and 2s1/2 in the ab- 
sence of multi-step process. Although other residual cen- 
tral interactions may also mix in such states, the mixed wave 
function from the tensor interactions is expected to have high- 
momentum components that are more significant than the 
others. 


We will measure the differential cross sections of (p,d) re- 
actions on '°Q. The advantage of this reaction lies in the se- 
lectivity of the momentum of the picked-up neutron. Com- 
bined with previous accomplished and acquired experimen- 
tal data [11-13], we would like to see the different behav- 
ior of the cross sections on the internal momentum of the 
picked-up neutron due to the effect of the tensor interactions. 
(p,d) reaction cross sections at 0° are measured for the tran- 
sition to the normal parity states, as well as to the abnor- 


mal parity state. Examples are the '°O(p, d)°O(gs +) and 
160(p, d)!SO(6.2 MeV 2.) as normal parity transitions and 


160 (p,d)O(5.2 MeV a”) as an abnormal parity transition. 
The effect of tensor force on the momentum distribution is 
expected to appear strongly in an abnormal parity transition. 
Therefore, the comparison of the energy dependence of the 
cross sections of those transitions can provide clear evidence 
of the tensor forces. Due to the large high momentum ampli- 
tude resulting from the tensor forces, the ratio of the cross sec- 
tion of abnormal parity state to the normal parity state would 
increase. 


Previous experiments were performed in RCNP WS- 
beamline at Osaka university [13]. We already observed 
a marked enhancement in the ratio for the a and/or a 
state(s) in 1O. However, because of the energy limit in RCNP 
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(400 MeV), the experiment has been performed at around 10 
degree to cover the required transferred momentum around 
2fm—!. To minimize other reaction mechanisms at 0 degree, 
we proposed an experiment of 400-1600 MeV/u !6O (p,d) 15O 
reactions at a deuteron angle 04 = 0° using the RIBLL2 sep- 
arator at HIRFL. 

To obtain the momentum distributions of nucleons in !6O at 
around 2fm~!, the proton beam from the CSRm synchrotron 
should be high quality and its energy range from 400 MeV 
to 1600 MeV in every 100 MeV step. The proton beam is 
injected into the target material, either a Carbon plate or ice 
target (H,O). We measured the (p,d) reaction for transitions 
to the ground state and low excited states, as well as the deep 
s-hole state. 


Il. SIMULATION PROCESS 
A. Proton beam and RIBLL2 setting 


The proton beam from the CSRm synchrotron should be 
high quality: the emittance in both directions should be less 
than 107 mm mrad (37 mm mrad would be the best), the mo- 
mentum broadening should be less than 0.02%, and the inten- 
sity of proton beam should be more than 10’ particles per pe- 
riod (10° if possible). To obtain the momentum distributions 
of nucleons in '?C and !6O at around 2 fm~!, the incident pro- 
ton energies will be changed from 400 MeV to 1600 MeV. In 
Table 1, we show the energy points in this experiment and the 
momentum transfer from the (p,d) reactions on 12C and 16O 
at 0°. 


TABLE 1. Momentun transfer for (p,d) reactions on 12C and 16O at 
0° 


Ep (MeV) Momentum transfer (fm7 +) 
2G 166 
400 1.50 1.54 
500 1.68 1.71 
600 1.83 1.86 
700 1.96 1.99 
800 2.07 2.11 
900 2.18 2:21 
1000 2.27 2.31 
1100 2.35 2.39 
1300 2.50 2.54 
1600 2.69 2.73 


The CSRm and RIBLL2 facilities are shown in Fig. 1. The 
proton beam, with a beam size around 1 mm, will hit the tar- 
get at FO, and the RIBLL2 will be run as a 0° spectrometer 
with a 75 mm (radius) drift. Since the neutron s-hole states in 
12C and '°O are expected to lie around 25 MeV and the typical 
width is 10 MeV, we plan to measure excitation energy up to 
about 40 MeV to cover the s-hole states. For this purpose, the 
slit at the first focal plane of RIBLL2 is set 6mm. At this 
configuration the momentum acceptance of RIBLL2 is +4%, 
which can cover the excitation energy up to 40 MeV, even at 
the lowest proton energy with 400 MeV. Given that the proton 
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will scatter from the edge of the slit, simulation of the multi- 
scattered proton has been done and only a small amount of the 
proton will come through the slit. Figure 2 demonstrates the 
transition of the proton beam and scattered deuteron before 
the target (upper-left), after the target (upper-right) at FO, af- 
ter the first quadrupole Q1, before the slit (bottom-left), and 
after the slit (bottom-right) at Fl. By properly setting the 
magnetic field of Q1, the primary proton beam and scattered 
deuteron will be separated clearly. 


B. Detector system 


In this experiment, we will monitor the intensity of the pro- 
ton beam using a detector system. The four silicon detec- 
tors are installed at FO in the target chamber and detect the 
light particles produced in the production target. Each sili- 
con detector has a dimension of 30mm x 30mm x 0.3mm. 
Figure 3 shows a schematic drawing of the beam monitor 
detectors. Through the counts of four silicon detectors, the 
position of the primary beam spot on the target can be moni- 
tored. The absolute calibration of the intensity can be done by 
the Faraday cup at the F1 focal plane or by the plastic scintil- 
lators at the F2 and F4 focal points. To maintain the vacuum 
in CSRm, a thin metal foil will be installed to isolate CSRm 
from RIBLL2. The energy loss and multi-scattered deuterons 
have been taken into account and simulated. 

For the detection and identification of the scattered 
deuterons from other particles, we use the time of flight (TOF) 
and energy loss (AE) method. Figure 4 presents the layout of 
detectors in this experiment. The TOF will be measured by 
two sets of plastic scintillation detectors (TOF 1 and TOF2) 
installed at F2 and F4 with a flight distance of 29.3m. The 
AE will be obtained from the TOF2 detectors. All the detec- 
tors in Fig. 4 will be placed in the atmosphere. The vacuum 
of RIBLL2 will be isolated by a 100 um thick Teijin DuPont 
film. Energy loss inside the film is negligible. This system 
allows us to identify protons and deuterons, and provides for 
separation of excited states in the final nucleus. 

The TOF 1 detectors are made of BC408 plastic scintillators 
with a dimension of 50mm x 50mm x 5mm. Two Hama- 
matsu H2431-50 photomultiplier tubes are used to read out 
signals from both ends. H2431-50 has good time perfor- 
mance: anode rise time is 0.7 ns and transit time is 16 ns. The 
energy losses in the TOF 1 detectors for the deuterons in all 
the experimental cases are about 1-3 MeV, which can pro- 
duce enough light for the detection of deuterons. The angular 
spread due to multiple scattering in TOF is estimated to be 
less than 2 mrad, which is negligible for our purpose. 

The TOF2 detectors consist of three BC408 plastic 
scintillation detector units. Each unit has a dimension 
300mm x 100mm x 10mm and is viewed by two Hama- 
matsu H2431-50 photomultiplier tubes. The TOF2 detec- 
tors have an active area of 300mm x 300mm. The large 
area is chosen to meet the beam spot requirement. The TOF 
vs. AE spectrum for the '°O (p,d)!°O reaction at an incident 
proton energy of 400 MeV has been simulated and shown in 
Fig. 5(a). In the simulation, the time resolution of the TOF 
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Fig. 1. (Color online) RIBLL2 and experimental setup. 
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Fig. 2. (color online) The transition of the proton beam and scat- 
tered deuteron before the target (upper-left), after the target (upper- 
right) at FO, before the slit (bottom-left), and after the slit (bottom- 
right) at F1. 


detector is 200 ps, the energy resolution of the AF detector is 
10%, the position resolution of the track detector is 0.4mm in 
horizontal and 0.5 mm in vertical, and the excitation energy 
of 15O is 5.18 MeV and 6.2 MeV. Here we select an ice tar- 
get (H,O) with a 0.46 g/cm? thickness as the oxygen target, 
which will be discussed later. 

Although the scattered deuterons can be identified from the 
TOF and AE method, it is hard to separate the excitation 
states of ''C and 150 due to the resolutions of the TOF and 
AE detectors. If we want to get good resolution of the ex- 
citation energy, track detectors (MWDCs) should be used to 
get precise momentum information of the deuterons. If we 
put the MWDC1, MWDC2, and MWDC3 at the focus points, 
according to the ion optics, the momentum P of the deuterons 
can be calculated by 


() 
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Fig. 3. (Color online) Schematic drawing of the beam monitor de- 
tectors. 
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Fig. 4. (Color online) Layout of detectors for deuterons. 


where Bp, Mẹ = —0.48 and D = 1.17m denote the mag- 
netic rigidity setting, the horizontal magnification, and the 
dispersion of RIBLL2, respectively [15]. 2x2 and x3 de- 
note the horizontal coordinate of the deuteron measured by 
MWDC2 and MWDC3, respectively. Figure 5(b) shows the 
simulation result of the momentum as a function of angles. 
The angular information of deuterons can be obtained from 
the position coordinates of MWDC3 and MWDC4. In the 
simulation, the position resolutions of the drift chambers are 
0.4mm in horizontal and 0.5 mm in vertical, respectively. It is 
clear that the excitation energy of '°O can be fixed. The focus 
of this method is that the MWDC1, MWDC2, and MWDC3 
must be placed at the focal points; otherwise, the resolutions 
of momentum spectra get worse. Therefore, the demand for 
the RIBLL2 beam tuning is quite high. 

In this experiment, the active areas of MWDC1, 
both of MWDC2 and MWDC3, as well as MWDC4 
are 150mm x 150mm, 300mm x 300mm and 
480mm x 320mm. Each MWDC has six detection planes. 
In total there are about 740 signals for MWDCs. For the 
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Fig. 5. (color online) The simulated results for the 160 (p,a) t50 reaction at an incident proton energy of 400 MeV. (a) Simulated TOF vs. 


AE spectrum, (b) simulated Momentum vs. angle spectrum. 


read-out of the anode signals and the time digitizing, REPIC 
RPA-130 preamplifiers and CAEN V1190A TDCs will be 
utilized. After the time digitizing, the data will be transferred 
to a computer through a VME crate. In order to obtain 
precise momentum and angular information, the position 
resolution of the drift chamber should be better than 0.4 mm 
in horizontal and 0.5 mm in vertical. The current simulation 
is based on simple calculations from a few hits with spatial 
coordinates (X, Y and U), and the set spatial resolution 
(0.4mm (X) and 0.5mm (Y )) is very conservative. A much 
better spatial resolution (< 250 um) shall be achievable with 
the CSR-RIBLL2 gas detectors, as demonstrated in previous 
published results [15, 16]. The feasibility of Fig. 5 is still 
on the safe side. Considering the angular resolution and the 
size of the beam spot at MWDC4, we set a 0.3m distance 
between MWDC3 and MWDC4. After taking into account 
the particle distribution at F4, both of the sizes of MWDC2 
and MWDC3 are selected to be 50mm x 50mm. 


C. Reaction target 


An ice target (H,O) with a 0.46 g/cm? thickness will be 
used in this experiment. The principle of choosing target 
thickness is to improve the yield ratio of deuteron under the 
condition of ensuring the resolution of the excitation energy. 
The effect of the energy loss and multiple scattering of the 
deuterons in the target at 400 MeV on spectral resolution have 
been simulated in Eq. (2). The angular distribution of mul- 
tiple Coulomb scattering is roughly a Gaussian distribution 
with a width given by 


Ee N aa h + 0.0381n (=)| [rad], 
Bep Xo 6 


here p, e, and z are the momentum, velocity, and charge 
number of the incident particles, respectively. x/Xo is the 
thickness of the scattering medium in radiation lengths. For 


the 460 mg/cm? 160 target, oy, is about 3 mrad (much less 
than emittance). On the other hand, the energy straggling 
caused by multiple scattering increases the momentum spread 
after the target, as a result, the energy resolution is deteri- 
orated. Figure 5(b) shows the momentum of the produced 
deuterons as a function of the emission angle. It can be seen 
that the first two excited states of !°O can be separated well 
at this thickness. 
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Fig. 6. The F1 simulated horizontal position spectrum of deuterons 
and scattering protons at a proton energy of 400 MeV. 


D. Background consideration 


The most likely dominant background in the (p,d) reactions 
is the elastically scattered protons. However, the momentum 
of the elastically scattered protons is different from that of 
the deuterons by about 20% at all incident proton energies. 
Figure 6 shows the simulation result of the horizontal position 
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spectrum at F1 of deuterons and scattered protons around 0° 
at a proton energy of 400 MeV. Since the momentum accep- 
tance of RIBLL2 is set as +4%, we expect to obtain deuteron 
spectra free from the proton background. 


Ill. RESULTS AND DISCUSSION 


At the proposed scattering angles of the deuterons, the cor- 
responding angles in the CM frame of the d +!4C or d +4°O 
system are smaller than 15°, i.e., larger than 165° in [17]. 
At this angle, the cross sections for the (p,d) reaction are ex- 
pected to be in the region of 20 ub/sr to about 100 ub/sr. As- 
suming the cross section to be ten times smaller than 20 ub/sr, 
and taking into account the beam intensity (J = 107 parti- 
cles per period), the target thickness (g/cm7), the transfer ef- 
ficiency of the beam line (€), and the angular acceptance of 
RIBLL2 (Q = 2.5 msr), the yields of measured deuterons are 
estimated as follows 


_TxKaxQxtx Na xe 


Y 
Mr 


(3) 


where Na is Avogadro number (6.02 x 1073) and Mp is 
molecular weight. Before further beam-optic simulation, and 
based on RIBLL2 beam-line design, we first assumed the 
transfer efficiency (€) to be 100%, and estimated the time re- 
quired to obtain about 2000 counts of deuteron events for each 
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energy. The estimated yields and the corresponding running 
time requirements are summarized in Table 2. 


TABLE 2. Estimated yield rates and time requirements 


Target Target thickness (g/cm?) Yield (period ty Time (h) 
C 0.45 1.12 100 
Ice(H,O) 0.46 0.77 140 


IV. SUMMARY 


In order to obtain good energy and angular resolutions, 
as well as a high production yield for the 400-1600 MeV/u 
160(p,d)!5O and '*C(p,d)!'C reactions at a deuteron angle 
64 = 0°, we used Monte Carlo method to simulate the kine- 
matic process based on the experimental conditions in the 
HIRFL-RIBLL2 separator. After taking into account target 
thickness, beam intensity, detector resolution, and so on, the 
TOF and AF detector settings were designed and the position 
and number of MWDCs were simulated. Total beam time and 
expected cross sections were also estimated. According to 
software simulation, some hardware conditions, such as beam 
monitor, high voltage optimization, and background product, 
need to be considered to achieve the expected results. This 
method will be continuously explored and developed from the 
experiment. The effects of tensor force can be systematically 
studied. 
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